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Particle number fluctuations are studied in the microcanonical ensemble. For the Boltzmann 
statistics we deduce exact analytical formulae for the microcanonical partition functions in the 
case of non-interacting massless neutral particles and charged particles with zero net charge. The 
particle number fluctuations are calculated and we find that in the microcanonical ensemble they are 
suppressed in comparison to the fluctuations in the canonical and grand canonical ensembles. This 
remains valid in the thermodynamic limit too, so that the well-known equivalence of all statistical 
ensembles refers to average quantities, but does not apply to fluctuations. In the thermodynamic 
limit we are able to calculate the particle number fluctuations in the system of massive bosons and 
fermions when the exact conservation laws of both the energy and charge are taken into account. 
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I. INTRODUCTION 

The statistical hadron gas model (see e.g. Ref. Q and recent review 0) appears to be rather successful in describing 
the data of nucleus- nucleus (A+A) collisions for particle multiplicities in a wide range of the collision energies. Usually 
one considers a thermal system created in A+A collision in the grand canonical ensemble (GCE). There are, however, 
situations when the canonical ensemble (CE) or even microcanonical ensemble (MCE) Q with explicit treatment 
of charge conservations or both charge and energy conservations are required. This happens, for example, when the 
statistical model is applied to elementary pp, pp, e + e~ collisions. Different statistical ensembles are not equivalent 
for small systems created in these collisions. 

In A+A collisions one prefers to use the GCE because it is the most convenient one from the technical point 
of view and due to the fact that both the CE and MCE are equivalent to the GCE in the thermodynamic limit 
(i.e. when the size of the system tends to infinity). However, the thermodynamic equivalence of ensembles means 
only that the average values of different physical quantities calculated in different ensembles are equal to each other 
in the thermodynamic limit. On the other hand, the analysis of fluctuations is also an important tool to study a 
physical system. Event-by-event analysis of A+A collisions (see e.g. Ref. |5() can reveal new physical effects not seen 
in observables averaged over a large sample of events. An essential part of the total fluctuations measured on the 
event-by-event basis is expected to be the thermal fluctuations. It was demonstrated for the first time in Ref. [f| that 
particle number fluctuations are different in the CE and GCE even in the thermodynamic limit. In the present paper 
we extend the results of Ref. [6j and make the analytical calculations of the particle number fluctuations in the MCE. 
In the textbooks of statistical mechanics the particle number fluctuations in the CE and MCE were not considered 
since the discussion was usually limited to the non-relativistic cases and the number of particles are assumed to be 
fixed in these ensembles. In the relativistic situation one can only fix the conserved charges (in the CE), or both 
the energy and conserved charges (in the MCE), while the particle numbers still fluctuate both in the CE and MCE. 
Results of both Ref. Q and the present paper demonstrate that the particle number fluctuations are different in 
various ensembles even in the thermodynamic limit. 

The paper is organized in the following way. First, we deduce and study the exact analytical expressions in the MCE 
for massless neutral particles with Boltzmann statistics: the partition function and average number of particles in 
Sec. II and particle number fluctuations in Sec. III. The extension of these results for the system of charged particles 
with zero net charge is considered in Sec. IV. In Sec. V we use the method of microscopic correlator proposed in 
Ref. [3- This gives us a possibility to study the MCE fluctuations in the thermodynamic limit for a much more 
general situation. It includes the effects of quantum statistics and non-zero particle mass. We consider the system 
of neutral particles and then extend the formulation to charged particles too. Both the energy and charge exact 
conservation laws are imposed. We summarize our consideration in Sec. VI. 
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II. THE MICROCANONICAL PARTITION FUNCTION AND AVERAGE NUMBER OF PARTICLES 

In order to calculate analytically the microcanonical partition function 1 we start in Sees. II-III with the system of 
non-interacting massless neutral particles and neglect the effects of quantum statistics (the extensions will be treated 
in Sees. IV- V). The microcanonical partition function for one massless particle with energy E in the volume V can 
be easily calculated: 

W 1 (E,V) = ^Jd 3 pS(E-\p\) = ^-^fdpp'SiE-p) = ^E*. (1) 

Here g is the degeneracy factor. In the case of two massless particles: 

1 gV 



W 2 (E,V) 



2V2.V Jo dqq2{E - q)2 =w{^) < 2 > 

The factor 1/2 appears because particles are identical. It can be proven that the ./V-particle microcanonical partition 
function Wn{E, V) has the form (see Appendix A): 

1 x N 

WN ^ = E (w-iy.Nr (3) 

where 

x = gVE 3 /n 2 . (4) 

The total partition function in the MCE is: 

oo oo jy oo n 

W(E,V) s £ w K (E,V) - s Y. j^W. = f £ (STWTT)! < 5 » 

N=l N=l y ' n=0 v ' v ; 



x ^ 3x n x ^ (x/27) r 



£^ (3n + 3)!n! 2£ ^ (*) (|) (2) 
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where we use (see e.g. Ref. and Appendix B) the Pochhammer symbol [a) n 1B2|I and the generalized hypergeometric 
function (GHF) p F q (ai, a p ; &i, b q ; z) (|Bl|l . 

In the GCE the independent variables are the volume V and temperature T. For neutral massless particles with 
Boltzmann statistics one easily finds the average number of particles and average energy: 

= <£>,„,. =MZ1. (6) 

7T 71" 

We want to compare the results of the MCE and GCE at equal volumes V and energies (E) g , c , e , = E. From Eq. © 
it follows: 

(iV> g . c ., = = (-) . (7) 



1 We define the microcanonical ensemble as the statistical system with fixed energy. An exact conservation of momentum, angular 
momentum and parity are neglected here. The MCE with conserving momentum was considered in Ref. 
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FIG. 1: The ratio of average particle number (N) m . c . e . © in the MCE to that N Q in the GCE. 



The average number of particles in the MCE equals to: 
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From Eq. (JSJ, using the asymptotic behavior of the GHF (|B3|) . one finds the asymptotic expansion of (N) m . c , e . © 
at N — > oo: 

Wm.ce. ^N(l + -±=+ — ^ + ..) . (9) 

V 8 N H52 N J 

The MCE and GCE are thermodynamically equivalent, i.e. (N) m . c . e . — N, in the thermodynamic limit VE 3 — > oo. 
The dependence of the ratio (N) m , c . e ./N on iV is shown in Fig.Q] 

As seen from Fig. ^the ratio (N) m , c . e ,/N goes to 1 at AT ^> 1. As the smallest number of particles in the MCE 
equals to 1, one finds that (N) m . c , e , ~ 1 when the volume V and the energy E of the system become small, i.e. 
if F_E 3 <C 1 ■ On the other hand, in the GCE one finds (N) g . c , e . cx (Vi? 3 )* ^ 1 in this limit of small systems. 
Therefore, aX N <\ the ratio (N) m c e /N is larger than 1 and it increases monotonously when N — > 0. 



III. PARTICLE NUMBER FLUCTUATIONS 



To study the particle number fluctuations in the MCE we calculate 
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The asymptotic expansion of (N 2 
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1(10)1 can be found using Eq. (IB3II 
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A measure of the fluctuations, the scaled variance w, is defined as usual, 

_ (iV 2 )„,c.e. ~ (N) 2 m . c . e . 



and it is plotted in Fig. [21 

The asymptotic expansion of ui m . c .e. (|12|l has the form: 
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FIG. 2: The scaled variance uj m .c.e. tHJl in the MCE. 
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FIG. 3: The particle number distributions in the MCE JT2J and GCE JT^ for N = 0.5 
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FIG. 4: The same as in Fig.0 but for TV = 10 
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In Figs. |3 and 0] we present the particle number distributions in the MCE and GCE for small (N = 0.5) and large 
(N = 10) systems. 
In the MCE the particle number distribution equals 

P (E V N) Wn(E > V) I X - N=l 2 (14) 

Wl 1 ' ' ~ W(E,V) ~ E W(E,V) (37V- 1)!JV! ' [ ' 

and it has the Poisson form in the GCE 

— N 

P g .c e . (N, N) = exp (- N) — , N = 0, 1, 2, . . . (15) 
Note that Poisson distribution P g . c . e . (N, N) (fl5|l results in 

W - N 2 

pg.c.e. = = = 1 ■ (16) 

It is seen from Eq. 1(13(1 that in the thermodynamic limit V — > oo the MCE scaled variance equals to uj m .c.e. = 1/4, 
and it remains quarter as large as the scaled variance u> g . c .e. GEt in the GCE. We conclude that the particle number 
distributions are different in the MCE and GCE. The MCE particle number distribution l|14|l is narrower than that 
of the GCE lfT31) in both small (N < 1) and large (N > 1) systems. At N < 1 the probability distributions 
Pm.c.e. (E, V, N) 1(141) and P g . c . e .(N, N) 114fl both have their maximums at the smallest values of N. Then the crucial 
difference between two ensembles follows from the fact that the minimal allowed value of JV is N = in the GCE, 
but N= 1 in the MCE (see Fig. 01. 

In the thermodynamic limit the average number of particles goes to infinity and the main contribution to the 
microcanonical partition function W(E, V) (JSJ comes from the states with large number of particles N 3> 1. In this 
limit the MCE particle number distribution can be simplified in the following way. Using the Stirling formula for 
factorials one finds: 

W N (E,V) = ^ {3N f iy m ^ex P [ /( A0], (17) 

where f(N) ~ TV log (Jbj — 4/V(log N — 1) + \ log 3 — log(27r) . We expand the function f(N) in Taylor series near 
the point of its maximum TV. For N ^> 1 and \N — N\ -C N we find: 

f(N) ~ f(N) + f"(N) \ (N-N) 2 = f(N) - 2 ( N ~ N ^ , (18) 

where N = (x/27) 1 / 4 is fixed by the condition f'(N) = 0, and it coincides with the result of Eq. Q. From Eq. 118(1 
it follows that the particle number distribution in the MCE can be approximated as: 



P m .c.e.(E,V,N) cx exp 



2 (N-N)' 



(19) 



For the Gauss distribution Pq(N) cx exp \—(N — 7V) 2 /2cr 2 ] the variance is easily calculated at N — * oo and it equals 
to (iV 2 ) - (N) 2 = a 2 . Therefore, from Eq. (T£J it follows (7V 2 ) m . c . e . - (N) 2 m c e = iV/4, and the MCE scaled variance 
is LOm.c.e. = 1/4. The Poisson distribution 1(15(1 for N ^> 1 can be also approximated by the Gauss distribution, but 
it equals to P g . c . e . (N, N) cx exp [-(N - iV) 2 /2lV] , and this leads to ui g . c . e . = 1. Therefore, at N > 1 both the MCE 
and GCE particle number distributions can be approximated by the Gauss distributions with the same average value 
N, but with different widths: u 2 n c e = N/4 and a 2 ce = N. A consequence of this is that the MCE scaled variance 
is a quarter the size of that in the GCE for classical massless neutral particles in the thermodynamic limit. 



IV. THE MCE FOR MASSLESS CHARGED PARTICLES 



The microcanonical partition function discussed in the previous sections can be generalized for the system of charged 
particles. If the system net charge Q equals zero the number of positively charged N + and negatively charged iV_ 



particles are equal in each microscopic configuration. The total MCE partition function is (a = gV/ir 2 )- 

OG OO n QQ /»00 

W(E,V,Q = Q) = dE + dE_W N+ {E+,V)W N _{E^,V) 

Ar + =l JV_=1 ^° ^° 

oo /*°° 

x5(N + -N_)5[E-{E + +E_)} = V / dE+ W N+ (E+, V) W N+ (E - E+, V) 

N+=l Jo 

- 1 m N+ , /v^-^-^r-. 

The last integral in Eq. 112011 can be easily evaluated using the Euler Beta-function: 

? w + -i (3iV + — l)! 2 
{6N+ - 1)! ' 



f dE+ E 3 + N+ -\E - E + ) 3N +~ 1 = E 6 ^- 1 B{3N+, 3N+) = E 6 
Jo 



Finally one finds: 



W(E,V,Q = 0) = J£ (6JV + -1)!JV+P = E^ (6JV+- l)!JV + !=» 



¥^ (6n + 5)!(n+l)! 2 ~ 120£ 7 I' 6' 3' 2' 3' 6 ' ' ' V216/ 



n=0 



3n + 5)!(n + l)! 



Similarly to Eqs. ((SJ ll(J|) after some calculations one obtains: 

p (, I 7 4 3 5 11 9. / g 

0^7 gj 31 2' 3' 6 ' ' V2167 ^ 

PL (■ 1 4 3 5 11 9. f' z V 
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±/m.c.e. 



and 



0^71) J-;- 1 -; 6 j 3 j 2 ' 3' 6' 1 216 J 



(-^±)m.c.e. - , 

P /. 7 4 3 5 U 9 9. / a y 
7 I' 6' 3' 2' 3' 6 ' ' ' \216/ 

In the GCE of Boltzmann massless charged particles with Q = (N + ) g c e ^ — (iV_) 9 . c .e. = one finds: 

gVT 3 QgVT A 

\N±)g.c.e. = — , (E)g.c.e. = 2 ' 

7T 7T 

The results of the MCE and GCE are again compared at equal volumes V and energies (£ , ) s . c .e. = E 
Eq. it follows: 

(iV ± ), c .e. = iV ± = (_) . 

The asymptotic expansions for (N±) m . c . e . (|23[) and (A f |), I1 . c . e . <|24|) at N± — ► 00 are found using Eq. (|B3fl : 

49 49 \ 



and 



(iV±)m. C . e . — N± 1 + _, 

V 23047V ± 9216iV± 



l ,T2^ ^ L 1 49 49 
(iVl) m .c.e. * N ± 1 + _ + ^ + =qr 

V 8iv ± 11527V ± 6144iV ± 



The behavior of the ratio (N±) m , c , e J N ± is shown in Fig. [SJ 
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FIG. 6: The scaled variance u>m. c .e. iHH f° r positively (negatively) charged particles in the MCE at zero net charge Q = 
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FIG. 7: The particle number distributions of positively (negatively) charged particles in the MCE 1311 and in the GCE (the 
Poisson distribution) at zero net charge Q = and N± — 0.5. 



The behavior of the scaled variance, 



rn.c.e. 



is shown in Fig. The asymptotic expansion of i 



{N±) m .c.e. 

at N± — ► oo has the form: 



1 



49 



while 



= 1 



11527V, 



(29) 



(30) 



. , , , — , , the same as for neutral particles, and uo^.e. = 1/2 
= 1/8 in the thermodynamic limit, and this is by a factor of 1/i 



6] . From Eq. H3U|) and Fig. HJ1 one sees that 
smaller than the scaled variance in the GCE 



and by a factor of 1/4 than in the CE. Therefore, for the system of massless particles with Boltzmann statistics the 
exact energy conservation leads to the MCE suppression of the particle number scaled variance in the thermodynamic 
limit by a factor of 1/4, and the exact charge conservation makes an additional suppression by a factor of 1/2. 

From Eq. (|2*2)l we find the positive (negative) particle number distribution in the MCE at net charge Q equal to 
zero: 



P m . c . e .(E,V,N±,Q = 0) = 



1 



r 2N ± 



E W(E, V,Q = 0) (6N±-1)\N±\ 



(31) 



In Figs.[7|and[S]we present the particle number distributions P m c e ,(E, V, N±, Q — 0) jSTJ for N± = 0.5 and N± = 10, 
respectively, and compare them with the GCE Poisson distributions which remain the same as for the case of neutral 
particles. Similar to Eqs. (|17I19|) we find the Gauss approximation for the MCE particle number distribution in the 
thermodynamic limit: 



P m . c . e .(E,V,N ± ,Q = 0) cx exp 



4 (N± 



lv±) 2 



(32) 



10 



0.4 



0.3- 



0.2- 



0.1- 











Q=0 



1 I 1 I 1 

■ 


m.c.e. 




• g.c.e. 


■ 

■ 

' i 




^ •■ •• f \ 

n -' ■ "0. 



10 

N 



15 



20 



+ 



FIG. 8: The same as in Fig. [7| but for N± = 10 



from which it evidently follows (N±), 



(N±) 2 m .c.e. = N±/8 and w± te . - 1/8. 



V. PARTICLE NUMBER FLUCTUATIONS FOR BOSONS AND FERMIONS 

To study the effects of quantum statistics in the MCE we use the technique proposed in Ref. Q. This method 
allows us to calculate the particle number fluctuations in the systems with the exact conservation laws imposed in the 
thermodynamic limit V — > oo. We reproduce the results of the previous sections for massless particles with Boltzmann 
statistics and study the MCE particle number fluctuations in a general case of the system of massive charged particles 
with Bose and/or Fermi statistics taken into account. 

Let us start with a system of neutral bosons or fermions 2 , then we extend our formulation to a system of charged 
particles. The system of neutral non-interacting identical Bose or Fermi particles can be characterized by the occu- 
pation numbers n p of single quantum states labeled by momenta p. The occupation numbers run over n p = 0, 1 for 
the fermions and n p = 0, 1, 2, . . . for the bosons. The GCE average values and fluctuations of n p equal to 0: 



•■■plg.c.e. 



((An p ) 2 ) 9 . 



g.c.e. 



exp(e p /T) - 7 ' 

— ( n p)g.c.e. — ( n p)g.c.e. 



{n p )g.c.e. (1 + l{np)g. c .e.) 



where 7 — +1 and 7 = 
Note that 7 = in Eqs. 



= vl 



(33) 
(34) 



•1 for Bose and Fermi statistics, respectively, e p = 



and to is the particle mass. 



corresponds to the Boltzmann approximation which is valid if (n p ) g . c . e . <C 1 for all p. 



2 There are examples of neutral bosons, like photon, n°, p°, etc., which are identical to their antiparticles. For fermions such a consideration 
has only illustrative purposes. One always needs to introduce some kind of charge to distinguish fermions from their antiparticles. 
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Expressions I|33I34[1 are microscopic in the sense that they describe the average values and fluctuations of single 
modes with momentum p. However, the fluctuations of macroscopic quantities of the system can be determined 
through the fluctuations of these single modes. To be more precise, we will demonstrate that the particle number 



fluctuations can be written in terms of the microscopic correlator (An p Ank) 
correlator can be presented as: 



g.c.e., where An p = n r 



Ig.c.e. 



This 



(An p An k )t 



v p 5 p k 



(35) 



due to the fact that the GCE fluctuations of the occupation numbers for different quantum states p ^ k are statistically 



independent. The variance ((AN) 2 ] 
to: 



g.c.e. 



(N 2 



g.c.e. 



(-W)g.c.e. of the total number of particles, N = n p , equals 



((AN) 2 ), 



E 

p,k 



'An p An k ) c 



(36) 



We have assumed above that the quantum p-levels are non-degenerate. In fact each this level should be further 
specified by the the projection of a particle spin j. Thus, each p-level splits into g = 2j + 1 sub-levels. It will be 
assumed that the p-summation includes all sub-levels too. This does not change the above formulation because of 
statistical independence of these quantum sub-levels. The degeneracy factor g enters explicitly when one substitutes, 
in the thermodynamic limit, the summation over discrete levels by the integration: 



£ 



2n 2 



p 2 dp .. 



(37) 



The scaled variance i 



'g.c.e 



in the thermodynamic limit V — > oo reads: 



_ ((AN 2 )), 



y v 2 



ig.c.e. 



(N)g. 



g.c.e. 



lg.ee. 



<■plg.ee. 



j™p 2 dpv 2 p 
Jo° C 'p 2 dp (n p ) g ., 



(38) 



The formula for the microscopic correlator will be modified if we impose the exact conservation laws in our equili- 
brated system. We introduce the equilibrium probability distribution W(n p ) of the different sets {n p } of the occupation 
numbers. In the GCE each n p fluctuates independently according approximately to the Gauss distribution law for 
An p with mean square deviation v 2 : 



W(n p ) cx J] exp 



(Anp) 2 



(39) 



To justify Eq. (|39|l one can consider (see Ref. 0) the sum of n p in small momentum volume (Ap) 3 with the center at p. 
At fixed (Ap) 3 and V — > oo the average number of particles inside (Ap) 3 becomes large. Each particle configuration 



inside (Ap) 3 



consists of (ApY 

,2 



gV/(2Tr) 3 » 1 statistically independent terms, each with average value 



'P/9-c.e. 

(|33|l and variance v 2 l|34l) . From the central limit theorem it then follows that the probability distribution for the 
fluctuations inside (Ap) 3 should be Gaussian. In fact, we always convolve n p with some smooth function of p, so 
instead of writing the Gaussian distribution for the sum of n p in (Ap) 3 we can use it directly for n p . 

Now we want to impose the exact conservation laws. The conserved quantity A (the energy and/or conserved 
charge) can be written in the form A = J2 p a(p) n p . An exact conservation law means the restriction on the sets {n p } 
of the occupation numbers: only those sets which satisfy the condition AA = y p a(p)An p = can be realized. Let 
us consider an exact energy conservation. Then A — > E (i.e. a(p) — > e p ) and the distribution (|39|l will be modified 
because of the exact energy conservation as: 



W(n p ) oc J^cxp 



p J V P ) p 



(An P Y 
2i 



exp 



, a + «A e p An p 



(40) 



where 6 (e p An p ) is the Dirac delta-function. It is convenient to generalize distribution (|4U|) using further the 
integration along imaginary axis in A-space. After completing squares one gets: 



W(n p ,X) oc Y[ 

exp 



(An p - Xv 2 e p ) 2 2 

+ 2 v p 6 p 



2v 2 



(41) 
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and the average values (i.e. the MCE averages) are now calculated as: 

= r^dXJZU.dn,... W(n p ,X) 

Using Eq. (|^|) one easily deduces 

({An p - v p Xe p )(An k - v 2 Ae fc ))m. c . e . = 8 pk v 2 p , (A 2 ) m . c . e . = - ^y^e^ , ((An p - v p Xe p )X) m .c.e. = . 

Therefore, one finds the MCE average for the microscopic correlator 

(An p An k ) m . c .e. = S pk - v p e p v\t k (A 2 ) + (An p X) v\t k + (An k X) v 2 e p (43) 

2 2 

r ,2 2 /\2\ r 2 ^P C P V k e k 

= pk + v p e p v k e k (A ) = d pk v p - — ^ 2 



By means of Eq. (|43[) one obtains: 



2 



((AN 2 )) m . c .e. S P , fc (An P An fc ) m . c . e , £ ^ 

Wm.c.e. = TTTT — 7 T — ^ , T _ V=W \ 5~~ 9 V 44 J 



(N) m . c . e . ^2 p {n p ) m .c.e. ^2p( n p) g.c.e. Y^p( n p) g.c.e. ^2 p v p e p 

j™p 2 dpv 2 p (f™p 2 dpv 2 e 



f™p 2 dp (n p ) g. c .e. f™P 2 dp (n p ) g . c . e . f™p 2 dp ^e 2 ' 

Comparing Eq. I|43[) and Eq. 1351) one notices two changes of the microscopic correlator due to the exact energy 
conservation. First, the MCE fluctuations of each mode is reduced, i.e. the value of ((An p ) 2 ) m . c .e. calculated from 
Eq. l|4*3|l at p = k is smaller than that of ((An p ) 2 ) g. c . e . = v 2 , given by Eq. l|3~4l in the GCE. Second, in the MCE the 
anticorrelations appear between different modes p ^ k (they are absent in the GCE). Both these changes result in the 
MCE suppression of the scaled variances (|44|l in a comparison with those in the GCE l|38l) . In fact, the first term in the 
r.h.s. of Eq. I|44|) equals to the GCE scaled variance w fl . c . e . I|38|l and the second negative term corresponds to the MCE 
suppression effects. Note also that according to Eq. 14411 the MCE fluctuations in the thermodynamic limit V — > oo 
can be presented in terms of the GCE quantities. The exact energy conservation should also lead to the differences 
between (n p ) m . c .e. and (n p ) g . c . e .. However, in the thermodynamic limit it follows: ^2 p (n p ) m . c .e. — ^2 p ( n p) g.c.e. ■ It 
can be also proven by the straightforward calculations that differences between (n p } rn . c .e. and (n p ) ff . c . e . lead to the 
corrections of the order of N to both (N 2 ) m . c .e. and (N)^ c e , but these corrections are equal to each other and they 
are cancelled out in the calculation of ((AN) ) m . c . e .. 

The GCE H38|) and MCE (|44|l scaled variances for different statistics are shown as functions of m/T in Fig. EI in 
the Boltzmann approximation (7 = 0) from Eq. I|38|) one finds: 

"sZt z = 1 , (45) 

which coincides with Eq. i|16|) used in Section III for m = 0. The Eq. (|45|l remains valid for all values of m/T. This is 
due to the fact that it follows from the Poisson particle number distribution P g B ° l ^ z (N, N) in the GCE, which is given 
by Eq. (|15l) at all values of m/T (only the average value of particle number TV decreases with increasing of m/T). 
From Eq. (|38|1 it follows that the effects of quantum statistics lead to the Bose enhancement, w^° s e e > 1 at 7 = 1, and 
the Fermi suppression, w^ 1 "™ < 1 at 7 — —1, of the particle number fluctuations. The strongest quantum statistic 
effects correspond to the m/T — > limit: 

^c.e e (™ = 0) = ~ 1.368, (46) 



9C(3) 



= 0) = ^ * 0.912 . (47) 



For the particle number fluctuations, as seen from Fig. EI the Bose enhancement ui Bose / LlJ Boltz > l and the Fermi 
suppression ui Fern ™ / \o Boltz < \ factors decrease monotonously with increasing of m/T, in both the GCE and the 
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FIG. 9: Three upper lines present the GCE scaled variances ui g . c .e. 1381 1 at different values of m/T, whereas three lower lines 
correspond to the MCE scaled variances u) m . c . B . 14411 . The dashed lines correspond to the Bose statistics, the dashed-dotted 
lines to the Fermi statistics and the solid lines to the Boltzmann approximation. 



MCE. These effects of quantum statistics in both ensembles become negligible at m/T 3> 1, as in this limit one finds 



( n p)c 



exp(— tp/T) <C 1, so that v 2 



for both the Bose and Fermi statistics. 



For the Boltzmann approximation all momentum integrals in Eq. (|44|l for 



p 2 dp exp = Tm 2 K 2 {- 



can be calculated analytically: 



T 



T 



2 7 / L P \ 

p dp e p exp ( -— ) = — 



«<® M?) 



> R 

P 2 dpel exp(-^) = - [K 5 (-) + K 3 (-) - 2K X (- 



Making use of the asymptotic behavior of the modified Hankel function K n (x) at x — > (i^o(^) 
K n {x) ~ \ T(n) (f )"" for n > 1) one gets in the m/T -> limit: 



<°i* e z (m = 0) 



1 

4 ' 



In x and 



(48) 



i.e. for classical massless particles the MCE scale variance is quarter as large as the corresponding scaled variance 
in the GCE. This result coincides, of course, with that of Eq. Ill.'ifl obtained in Section III from the MCE partition 
function of massless particles with Boltzmann statistics. For the case of Bose and Fermi statistics we obtain: 



1 (m = 0) = 



6£(3) 



9 £(3) 



135 £(3) 

2tt 4 
405 £(3) 

7tt 4 



0.535 , 
0.198 . 



(49) 
(50) 
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The Bose enhancement and Fermi suppression of the fluctuations exist in both the GCE and the MCE. However, the 
effects of quantum statistics for the particle number fluctuations are stronger in the MCE than those in the GCE . As 
it follows from Eqs. (|48H50(1 the Bose enhancement of the scaled variance in the MCE at m/T — > is approximately 
equal to 2.142, and the Fermi suppression 0.793. These numbers can be compared with 1.368 and 0.912, respectively, 
found from Eqs. gSHSJl in the GCE. 

As seen from Fig. ED the MCE suppression of the particle number fluctuations for massive particles is stronger than 
that for massless ones, and all MCE scaled variances u) m .c.e. decrease monotonously with increasing of m/T. From 
the asymptotic expansion at a; > 1 [j| 



K n (x) ~ \J ir/2x exp(— x) 



1 



^ (2x) k 

k>0 y ' 



(4n 2 -l)...[4n 2 -(2fc-l) 2 
k\2 2k 



(51) 



one finds the behavior of the scaled variance uJm^e — %{ m /T)~~ 2 < 1 at m/T 3> 1 for the Boltzmann approximation 
(to obtain this result one needs to keep terms up to k = 3 in Eq. J5TJ). The MCE scaled variances for the Bose 
and Fermi statistics have the same behavior, as the effects of quantum statistics are negligible at m/T ^> 1, so that 



, .Bose ^ . .Fermi 
^m.c.e. — ^m.c.e. 



0J 



Boltz 



For the system of charged particles and antiparticles in the GCE, similar to Eqs. (|33I34() . one has: 



((A< 



p Ig.c.e. 

±\2\ 



exp[(e p =F j") /T] 



Ig.c.e. 



( n p) 2 g.e.e. = (»p }g.c.e. (l + l{n p )g.e.e) 



,±2 



(52) 
(53) 



where /i is the chemical potential connected with the conserved charge Q. The microscopic correlator (I35f) can be 
then generalized as: 



(An"Ar 



k/gce. 



3pk &a{: 



(54) 



where a, (3 are + and (or) — . The average values of the energy E = J2 p a e p n p and charge 3 Q = ^Z pa q a n p are 
regulated in the GCE by the temperature T and the chemical potential fi, respectively. Similar to Eq. Q40[l the MCE 
distribution for the occupation numbers can be presented as: 



W(n p ) 



OC J^J GXp 



2vf 



\p,a / \p,a 



(55) 



OO 



dXj 



d\Q ]J[ exp 



^ff- + i\ E e p An% + i Xq q a An p 



where 5 (^J2 P a e p^ n pj an d ^ fSp a Q a An'^j correspond to the exact energy and charge conservations, respectively, 
in the MCE. After the straightforward calculations, similar to Eqs. 141H43|) . one gets an expression for the microscopic 
correlator in the MCE for the charged particles with the exact charge conservation law imposed: 



(An°An£) 



a2 /32 
V P V k 



*VE< e 



E<v 2 

p,oc 



(56) 



where 



w - E« 




E< 2£ p9° 

V P,ot 



(57) 



3 In what follows we assume g+ = 1 and q = —1, so that Q = jV+ — However, other values with q+ = — q can be also used. 
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Therefore, one finds: 

a _ ((AN*)) m . c . e . _ £ p , fc (An«A<) 



a) m.c.e. L^pV^p /m.c.e. 



(58) 



^2p{ n p) g-ce. ^2p( n p) g.c.e. ^2p,a V p 2e p Ep( n p) g.c.e. JJp.a V p 2( f 2 

The first term in the r.h.s. of Eq. (|58|l corresponds to the GCE scaled variance ojg Ce for the positive (a = 1) or 
negative (a — —1) particles. Note that in the Boltzmann approximation ojg t Pf ltz coincides with the scaled variance 
for the neutral particles oj^° lt e z l)45|l and, therefore, equals 1. The second and third terms in the r.h.s. of Eq. 1581) 
correspond to the MCE suppression of the fluctuations due to the exact conservations of energy (compare to Eq. (|44|) ~) 
and charge, respectively. 

In the case of zero net charge Q — (this means [i = and, therefore, (?i+) g . c . e . — (n~) g .c.e.) one finds for charged 
particles: 



E p v± 2 1 I> P ±2 1 (E P 



v ± 2 e 



U m.c.e.(Q — 0) — ± ± ± ±22 — LUm.c.e. , (59) 

l^p\ n V 1 g.c.e. z l^p\ n V 1 g.c.e. z L-ip\ n P I g.c.e. l^p V P e p Z 

where uj m .c.e. in Eq. I|59|) corresponds to the MCE scaled variance Q44[l for the neutral particles with the parameters 
m and g equal to those of the charged particles. Thus, we came to the conclusion that in the MCE at Q = 
an exact charge conservation leads to the scaled variances of (negative) positive particles which are by a factor of 
1/2 smaller than the corresponding MCE scaled variances for neutral particles. This result is an agreement with 
the CE suppression of the particle number fluctuations found in Ref. [||: the scaled variance in the CE at Q = 
in the thermodynamic limit equals to cj^ e = | e (note that at fi = the GCE scaled variances for charged 
and neutral particles are equal to each other). In the case of Boltzmann statistics one finds for massless particles 
^m B c°e tz ( m = 0) = 1/8 , which coincides, of course, with the result of Eq. (|30|) obtained in Sec. IV. 



VI. SUMMARY 



We have studied the particle number fluctuations in the MCE. First, in Sec. II we have obtained the partition 
function for the system of non-interacting massless neutral particles with Boltzmann statistics. This allows us to 
calculate the particle number fluctuations for this system in Sec. III. These MCE fluctuations are suppressed in 
comparison to those in the GCE. In the thermodynamic limit the MCE scaled variance of the multiplicity distribution 
equals 1/4 and this is a quarter the size of that in the GCE. As a second step we consider in Sec. IV the system of the 
charged massless particles with zero net charge and calculate its partition function. This leads to the scaled variance 
equal to 1/8 in the thermodynamic limit. Thus, an exact charge conservation makes an additional suppression of the 
scaled variance by a factor of 1/2. In Sec. V we use an approach proposed in Ref. 0- This method does not work 
accurately for finite (small) systems, but it does allow us to calculate the correct values of the scaled variances in the 
thermodynamic limit V — > oo for a much broader class of statistical systems - massive bosons and fermions - when 
both the energy and charge exact conservations are taken into account. The effects of quantum statistics and non-zero 
particle mass lead to the significant changes of the MCE particle number fluctuations, and they are studied in details 
in Sec. V. We have also reproduced the limiting V — » oo behavior of the MCE scaled variances for massless particles 
with Boltzmann statistics obtained in Sees. III-IV from the exact analytical expressions for the partition functions. 
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APPENDIX A 

The microcanonical partition function Wn(E,V) can be recursively calculated from Wn-i(E, V) : 

W N (E,V) = / # p W * J * p *-*.& / ^ W ^-ElP«l) (Al) 



fe=i 



T/ - Ar " 1 

X 

(2tt)3 

Now we shall prove by induction method that for arbitrary N > 1 the microcanonical partition function has the form 
of Eq. J3J). It can be checked directly that Eq. at N = 1 and N = 2 coincides with (JJJ and J5J, respectively. 
Now assume that it is correct for N — 1 , 

*WS,V) = 2 ":l TTi fS^'V"* , (A2) 



(3JV-4)!(iV- 1)! V2tt 2 
and substitute (|X2|) into Eq. IjAljl : 



Wn(e,v) = ^ (3jV _ 4)!(JV _ 1)! ^ ) (^-,r- 4 (A3) 



1 /fff 



j (3N-4)l(N-l)lJ dpp2 ( E ~rt 



2 I TTi ^\3W-4 



2V V2tt 2 7 (3AT-4)!(AT- 1)! ,/ 
The last integral can be easily evaluated after substitution p = E — q : 



E r E OZ?3JV- 

2 „3iV-4 Z - C/ 



/ dpp 2 (E~p) 3N - 4 = f dq(E-q) 
Jo Jo 

Eqs. HA3IA4H result in Eq. ©. 



(3iV-l)(3iV-2)(3JV-3) 



(A4) 



APPENDIX B 

The generalized hypergeometric function (GHF), also known as the Barnes extended hypergcomctric functionis 
defined by the following series, 

T-1 / , , , a {ai)k(a 2 )k ■ ■ ■ {a p ) k z k 
p Fq(a-L,a2, ■ ■ . ,a p ;bi,b 2 , . . . ,b q ; z) = > , , ,, ; r , , (Bl) 



where (a)k is the Pochhammer symbol: 



(a) fc = r( ° + fc) =o(o+l)...(a + fc-l) . (B2) 
r(a) 



The asymptotic behavior of GHF at z — > oo is given at p=0 by 



V(9 + l)(27r)« 



1 + -v* + ( + ( *+^) 3 + " ' ' ' s JJ ' ! ' 
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where 



and 



Cl = -ig 2+ q (j + 2 \, (B4) 
2 24(g+l)' V ' 

g 2 (g + 2) 2 + 768 n g(g + 2) [5g 2 (g + 2) 2 - 288g(g + 2) + 10944] 



2304( g +l) 2 414720(g+ If 



1 



D 2 = B 2 - ^Bf, (B7) 



_3_ 

(gTl)"^ ' (g + 1) 2 

4 „ „ 6 „,„ 3 



D 3 = B 3 - 7 —^B 1 B 2 + j- ■ B%, (B8) 



Di = Bi - (^TT) SljB ' 3 + (^TI)^ 2 - (gTi) 1 ^' (B9) 



B m = l + J2(bj) m . (BIO) 

4=1 
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